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ABSTRACT 

H 2 D + is a primary ion which dominates the gas-phase chemistry of cold dense gas. Therefore it 
is hailed as a unique tool in probing the earliest, prestellar phase of star formation. Observationally, 
its abundance and distribution is however just beginning to be understood in low-mass prestellar and 
cluster-forming cores. In high mass star forming regions, H2D + has been detected only in two cores, 
and its spatial distribution remains unknown. Here we present the first map of the ori/io-H 2 D + k~ 
= li,o -> 1.1,1 and N 2 H + 4-3 transition in the DR21 filament of Cygnus-X with the JCMT, and 
N 2 D + 3-2 and dust continuum with the SMA. We have discovered five very extended (< 34000 AU 
diameter) weak structures in H 2 D + in the vicinity of, but distinctly offset from embedded protostars. 
More surprisingly, the H 2 D + peak is not associated with either a dust continuum or N 2 D + peak. 
We have therefore uncovered extended massive cold dense gas that was undetected with previous 
molecular line and dust continuum surveys of the region. This work also shows that our picture of 
the structure of cores is too simplistic for cluster forming cores and needs to be refined: neither dust 
continuum with existing capabilities, nor emission in tracers like N 2 D + can provide a complete census 
of the total prestellar gas in such regions. Sensitive H 2 D + mapping of the entire DR21 filament is likely 
to discover more of such cold quiescent gas reservoirs in an otherwise active high mass star-forming 
region. 

Subject headings: ISM: abundances — ISM: clouds — ISM: molecules — ISM: structure — radio lines: 
ISM — stars: formation - ISM: individual (Cygnus-X) 



1. INTRODUCTION 

In their evolution toward the formation of a young stel- 
lar object, low-mass cloud cores change their density dis- 
tribution until the so-called "pivotal" state is reached 
(e.g. Shu et al. 1987). This stage is considered the start- 
ing point for protostellar accretion. Thus, it is crucial 
to find and study objects close to the "pivotal" stage to 
unveil the initial conditions in the process of star for- 
mation. We believe to have found this class of objects: 
they are the so-called prestellar cores, characterised by 
central densities of ~ 10 6 cm -3 , centrally concentrated 
density profiles, and with clear signs of contraction in the 
central few thousand AU (Caselli et al. 2002a, b; Crapsi 
et al. 2004, 2005). They stand out in having large degrees 
of CO freeze-out and deuterium fractionation (Bacmann 
et al. 2003). In particular, they show strong emission 
in the Jk+.k- = li.o I1.1 line of or£/io-H 2 D + in the 
case of LI 544, the prototypical low-mass prestellar core 
(Caselli et al. 2003). In this object, the H 2 D + line has 
been mapped by Vastel et al. (2006) and found to be ex- 
tended over a region of ^5000 AU in projected radius. A 
H 2 D + survey was conducted in low-mass cores (Caselli 
et al. 2008) and revealed a decrease in fractional abun- 
dance of ort/io-H 2 D + in protostellar cores as compared to 
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starless cores. 

To reproduce the past observations towards low-mass 
cores, which imply abundances of H 2 D + close to 10~ 9 
w.r.t. H 2 , one needs to have species such as CO (which 
destroy H 2 D + ) almost completely frozen onto dust grains 
(e.g. Walmsley et al. 2004; Roberts et al. 2004). This can 
be achieved only if the dust temperature is low (Td us t < 
20 K) and the density is large (n[H 2 ] > 10 5 cm" 3 ). At 
these conditions, all species quickly condense out (with 
time scales ~10 9 yr (n[H 2 ]/cm~ 3 )~ 1 , leaving mainly He, 
H 2 and its daughter species (in particular H^" and its 
deuterated forms, including H 2 D + ) in the gas phase. 
Therefore, H 2 D + observations are a unique tool to trace 
and study the chemistry and kinematics of the high den- 
sity nucleus of prestellar cores, the future stellar cradle. 

That being said about low-mass star formation, the 
early phase chemistry in massive star forming regions 
is still poorly understood, mainly because the study of 
these regions have so far focused on molecular mate- 
rial associated with high mass protostellar objects (e.g. 
Beuther et al. 2002), which are actively altering gas and 
dust properties. Based on our wide-field (10 x 10 arcmin) 
census of clouds in 32 massive star-forming regions, har- 
boring UCHn regions, we discovered a new sample of 
massive pre-protocluster cores. Our program is known 
as SCAMPS (the SCUBA Massive Pre/Protocluster core 
Survey; Thompson et al. 2005). Many of the sources 
are seen as dark patches in infrared images of the re- 
gion. As a result these clumps must have dust tempera- 
tures below 30 K (as evidenced by mid-IR upper limits) 
and have masses of a few 100 to 1000 Mq. These mm- 
peak/mid-infrared absorption sources are very promising 
candidates for being high mass pre-protocluster cores. 
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We started investigating these pre-protocluster cores 
in more detail by observing them in NH3 with the Effcls- 
berg 100 m telescope, as well as in deuterated NH 3 with 
the IRAM 30m telescope(Pillai et al. 2007). On aver- 
age, the derived temperatures are around 15 K compared 
with average temperatures toward the previously stud- 
ied more evolved High Mass Protostellar Objects (HM- 
POs, Sridharan et al. 2002) of 22 K. Finally, depiction 
and very large deuterium fractionation, [NH2D]/[NH3] 
- 0.2, are found (Pillai ct al. 2007, Pillai et al. 2011), 
an enhancement of « 10 4 over the cosmic D/H ratio. 
With gas masses of the order of a few hundred M , 
high dcuteration, depletion, and cold temperatures, these 
sources feature many characteristics expected for massive 
pre-protocluster cores. Similar large values of the deu- 
terium fraction have been found toward quiescent cores 
embedded in Infrared Dark Clouds (IRDCs; Chen et al. 
2010, Fontani et al. 2011). However see also other re- 
cent deuteration studies in IRDCs (Miettinen et al. 2011, 
Sakai et al. 2012). 

With the JCMT and SMA we observed the ortho- 
H 2 D+(li i0 -li,i), N 2 D+ 3-2, and N 2 H+ 4-3 line towards 
three cores from our sample. Here we report detection 
of the ort/io-H2D + (l i.o-l 1.1) line and its large -scale 
emission in one of the three cores (the only core where 
ort/io-H 2 D + was detected) and subsequent mapping of 
the deuterated region in N 2 D + 3-2 line with SMA. 

2. SOURCE SELECTION 

The JCMT observations were initially targeted to- 
wards three objects from the SCAMPS survey with the 
brightest NH 2 D detections. In total 23.8 hours were 
spent on the project on the three sources. After few hours 
of integration on each of the three targets which did not 
result in a detection, we decided to perform deeper inte- 
gration on G81.74+0.59, since we detected the brightest 
NH 2 D line towards it in our previous IRAM 30m sur- 
vey (Pillai ct al. 2007). Here, we only discuss the ob- 
servations towards G81.74+0.59. The maps presented 
here correspond to clumps along the DR21 filament ~ 2' 
north of the well known DR21(OH) massive young stel- 
lar object. The DR21 filament itself is the most massive 
region in the entire complex. Our mapped region coin- 
cides with three bright 1.3 mm MAMBO dust continuum 
peaks as revealed in a large scale survey of Cygnus-X 
(Motte et al. 2007). For reference, this region encom- 
passes the massive dense core N53 of Motte ct al. survey 
(see also Bontemps ct al. 2010, Csengeri et al. 2011) and 
FIR 1/2/3 in infrared studies of this region (Chandler 
et al. 1993, Kumar et al. 2007). The embedded pop- 
ulation identified by Kumar et al. is shown as stars in 
Fig.l. 

3. SMA AND JCMT OBSERVATIONS 

JCMT N 2 H + and H 2 D + observations were carried out 
in service mode in several sessions starting August 2007 
to August 2008 with the HARP instrument on JCMT. 
HARP is the heterodyne multi-pixel array at JCMT 
equipped with 16 receptors separated by 30 arcsec with 
a foot print of 2 arcmin. The beam size of each receptor 
is « 14 ". Nyquist sampled maps were achieved using a 
jiggle pattern with a position-switch to move to the ref- 
erence. During most of our observations, 14 of the 16 
receptors were working. After finding evidence of emis- 



sion at the edge of our map from the initial observing 
session, we shifted the map center for the subsequent ses- 
sions. The pointing accuracy was on average < 3 ", and 
the median system temperature at the H 2 D + frequency 
was 579 K. The smurf software provided as a part of the 
Starlink package was used to regrid the data, mask the 
bad receptors, construct 3D data cubes, and co-add the 
scans. The data were then converted to fits and analysed 
in CLASS 6 . 

The Submillimeter Array (SMA 7 ) N 2 D+ and contin- 
uum observations were made with all eight antennas in 
the most compact configuration (sub-compact) in two 
tracks in August and September 2008. The observations 
were done in track sharing mode with several sources per 
track. For G81.74+0.59, we had four pointings separated 
at less than half a (primary) beam spacing to cover the 
same region as the JCMT HARP footprint at 230 and 
279 GHz respectively. The 230 GHz receiver was tuned to 
the CO (2-1) line in the spectral band sl3 of the USB for 
the N 2 D + observations. We used a non-uniform correla- 
tor configuration to have a higher spectral resolution of 
~ 0.2 kms -1 for the chunks covering the N 2 D + and CO 
line. The rest of the correlator was set to ~ 0.4 krns" 1 . 
The data were taken under very good weather conditions 
(< 1mm water vapor), and the typical system temper- 
atures were 100 K. The gain calibrator was MWC349a, 
~ 2 °away from the target position. Uranus and 3C454.3 
were used as the flux and bandpass calibrator, respec- 
tively. Flux calibration is based on SMA flux monitor- 
ing of the observed planet and is estimated to be accu- 
rate to within 20%. The primary beam at 230 GHz is ~ 
51". The data calibration was done in MIR, an IDL based 
calibration package, and then exported to MIRIAD and 
GILDAS for further data reduction which included con- 
tinuum subtraction, and imaging. The center of the fi- 
nal image of the mosaic, the 1 a rms noise level achieved 
in the continuum images, and the resulting synthesized 
beam for all observations is listed in Table 1. 

4. RESULTS 

4.1. H 2 D+ Clumps 

The detection of JCMT H 2 D + is at a modest signal-to- 
noise (S/N) with five confirmed S/N peaks (5, 6, 7, 7 and 
9 a) for a smoothed velocity resolution of ~ 0.2 kms -1 , 
therefore we obtained only source candidates. In Fig.l, 
the S/N map obtained (after baseline subtraction) is 
shown. To increase the S/N for the peak positions, we 
smoothed the original JCMT map using a symmetri- 
cal Gaussian PSF with a full-width at half-maximum 
(FWHM) of two pixels, evaluated over a length of 10 
pixels. 80% of the pixels in the inner 90 x 90 "of the 
map shown in Fig.l have rms noise between 0.013 and 
0.021 K. This resulted in our final map presented here 
with a FWHM beam of 19.7 ". The > 3cr detections 
in the S/N map were then confirmed by visualisation 
of spectra in the smoothed map, and subsequent S/N 
analysis on Gaussian fits to the spectra. Similarly, three 
other > 3<T detections (also seen in Fig.l for example 

6 http: / / www. iram . fr /IRAMFR/GILDAS / 

7 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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the contour between P4, PI and P3) were rejected by 
visualisation of spectra. These might well be true detec- 
tions, however requires validation with deeper observa- 
tions. Following these steps, we have confirmed five reli- 
able detections (labelled PI to P5) whose positional off- 
sets and Gaussian fit parameters (obtained with CLASS) 
is tabulated in Table 2. The N2H + (4-3) line parameters 
for data smoothed to the same spatial resolution is also 
given in the same table. The spectra towards the five 
positions for both lines are shown in Fig. 2. The spectra 
have been converted from antenna temperature to main 
brightness temperature for a main beam efficiency of 0.63 
(from the JCMT HARP web page). 

4.2. Spectral Line Features 

The ti2D + , N2H+ and N2D+ spectral line parameters 
towards the five H 2 D + peaks are reported in Table 2. We 
also provide the non-thermal line width for these species 
after taking away the contribution due to thermal broad- 
ening at 15 K. 

The H2D+ line towards even the brightest core (say 
PI) at offset (37 ", 44 ") is only modest at 0.18 K. 
Our brightest line is still fainter than the weakest H 2 D+ 
detection towards low-mass starless and prestellar cores 
reported so far (for a complete list see Caselli et al. 
2008). The LSR velocity towards the northern sources 
(PI, P2, P3) is systematically offset from those of the 
southern cores (P4, P5). We find a range in line widths 
of 0.8 < Au/kms" 1 < 1.6. The thermal line width 
(FWHM) at 15 K is 0.4 kms" 1 . Clearly, the line width 
has a significant non-thermal component as expected in 
such protocluster environments (also given in Table 2) 
. Though it is difficult to assess the line profiles at the 
given S/N, a double-peaked profile is observed towards 
P4. The line asymmetry might be due to different ve- 
locity components along the line of sight, absorption in 
the core envelope or infall motions (as argued by van der 
Tak et al. 2005 where a similar H2D+ profile is observed 
in L1544). Clearly, high spatial resolution observations 
are required to clarify this. 

The 4-3 line of N 2 H+ line is very bright (> 2 K) and 
broad (> 1.5 krns" 1 ) towards all the five offset posi- 
tions. This transition has over 40 hyperfine lines within 
< 5 kms" 1 . The N 2 D+ 3-2 transition also has hyperfine 
structure (hfs) spread over 7 kms^ 1 . However this is be- 
low our detection limit for all H 2 D+ peaks. Also, the S/N 
for N2D + spectra towards H2D + peaks is poor. There- 
fore, we refrain from any interpretation of the N2D + 
spectra. H2D+ is much lighter than N2D+ and N2H+ 
(factor of 7) and hence the thermal line broadening is 
significant. The thermal line width (FWHM) at 15 K 
is 0.41 kins" 1 for H2D+, much larger than 0.15 kms -1 
for N2H" 1 " and N2D+. However, in Cygnus-X, the non- 
thermal contribution is still dominant. As listed in Ta- 
ble 2, the intrinsic (i.e. taking the hfs structure into 
account) line widths of three species are different. Ex- 
cluding P4 (where the N 2 H+ and H 2 D+ line widths are 
equal), H2D+ line width is lower than that of N2H+ and 
might imply that deuterated molecules preferably trace 
quiescent gas. The mean difference in line width for the 
H 2 D+ - N 2 H+pair is 0.52(±0.48). i.e. the higher N 2 H+ 
line width with respect to H 2 D+ is only marginally sig- 
nificant. 



The N2H+ LSR velocity is consistently red shifted rel- 
ative to that of H 2 D+ and magnitude of this shift varies 
from 0.2 to 0.6 kms -1 . 

4.3. H 2 D+, N 2 D+ 7 N 2 H+ and Thermal Dust Continuum 

The dust continuum emission at high (1.3 mm SMA, 
~ 3.8 ") and coarser (1.3 mm MAMBO, ~ 11 ") angu- 
lar resolution is shown in panels A and B of Fig.l. The 
MAMBO clumps fragment further with significant struc- 
ture all along the filament. The subsequent panels in the 
figure show the gas distribution in three different trac- 
ers: H2D + , N2D" 1 " and N2H+. There is very good corre- 
lation between thermal dust continuum and N 2 H + emis- 
sion with the brightest N 2 H+ emission coinciding with 
the bright MAMBO peaks (panel D). However, there is 
no such one to one correlation of dust continuum with 
either N 2 D+ (panel C) or H 2 D+ emission (panels B, C 
and D). 

H 2 D+ is concentrated in five clumps distributed 
around the massive protocluster (panel B Fig.l). How- 
ever, no H2D + emission is seen towards any of the bright 
and dense dust cores seen with MAMBO/SMA. Con- 
versely, no compact continuum is detected for H2D + 
cores. The H2D+ distribution coincides with our SMA 
N 2 D + 3-2 observations, however their peaks do not. In- 
stead the N 2 D+ cores are distributed in further clumpy 
structures offset and partially surrounding the H2D + 
peaks. Furthermore, there is one H2D+ peak (PI) with 
no associated N2D + emission, and several N2D + cores 
with no corresponding H 2 D+ emission. As seen in the 
fourth panel of Fig.l, there is always extended N 2 H+ 
emission associated with H 2 D+ cores. 

4.4. ort/io-H 2 D + and N 2 H+ abundance 

The ort/io-H 2 D+ column density for the different cores 
have been estimated following Vastel et al. 2006 from 
the main beam brightness temperature, line width and 
assuming an excitation temperature (Tex) of 10 K. The 
H2D" 1 " line is optically thin in all cores (<C 1) and the 
assumed Tex also enters the optical depth estimation 
following equation 4 of Vastel et al. 2006. Therefore, we 
note that the derived H2D + column density crucially de- 
pends on the assumption of excitation temperature. A 
30% decrease in Tex lower than 10 K is poorly justi- 
fied. Pillai et al. (2007) estimate a ammonia rotational 
temperature of 18 K in the vicinity of P4 (within 40"). 
The presence of H 2 D + emission together with the lack of 
24^m emission towards the H 2 D + cores rules out much 
higher gas temperatures (3> 20 K). In their H2D" 1 " sam- 
ple of prestellar and protostellar cores, Caselli et al. find 
on average 3 K higher kinetic temperatures and higher 
Tex for their protostellar cores as their prestellar cores 
(8 K vs 11 K). Friesen et al. 2010 also find a higher (than 
prestellar cores) Tex of 12 K to be appropriate for their 
low-mass cluster-forming core in Ophiucus. 

MAMBO 1.3 mm data smoothed to the resolution of 
the final JCMT data have been used to determine H2 
column density for a 15 K kinetic temperature following 
Kauffmann et al. (2008). We have used a dust opacity of 
0.0102 cm 2 /g for dust grains with thin ice mantles and 
gas density n(H) = 10 6 cm _3 as in Ossenkopf & Henning 
(1994). The ort/io-H 2 D+ column density and the result- 
ing H 2 D+ abundances are reported in Table 3. The ratio 
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between column density towards cores studied here and 
those observed towards low-mass cores (summarised in 
Table 3 of Caselli et al. 2008) varies from 1 to 13. 

As mentioned before, the N 2 H + hyperfine structure 
contribute significantly to the broadening of the line. We 
fit the N 2 H+ hyperfines (METHOD HFS in CLASS) and 
FWHM so obtained is much lower than from a GAUSS 
fit (Table 2, Gerin et al. 2001). The optical depth, 
and therefore the excitation temperature which is usu- 
ally otherwise derived based on optical depth and line 
brightness is uncertain. Hence the choice of excitation 
temperature is somewhat arbitrary and is based on obser- 
vations in literature. For example, Fontani et al. (2008) 
derive excitation temperature close to 6 K value towards 
N 2 H+ cores in a high mass star forming region. As- 
suming optically thin condition, we have estimated the 
column density from the integrated intensity of N 2 H+ 
4 — 3 lines following Caselli et al. 2002c for an excitation 
temperature of 5 K. The total N 2 H + abundance is tabu- 
lated in Table 3. The uncertainties given in brackets are 
formal errors obtained by a Gaussian error propagation. 

5. DISCUSSION 

While some of the results discussed above are intu- 
itively easy to understand, many are not. In this section 
we attempt to explain these anomalies and thus charac- 
terise the physical conditions. 

The bright MAMBO sources are actively forming stars 
and clusters (Kumar et al. 2007), as also evidenced by 
their mid-infrared (mid-IR) emission. Higher temper- 
atures, and subsequent release of heavy neutrals (CO) 
from grain mantles destroy H 2 D + and therefore reduce 
the deuterium reservoir in such bright cores. For ex- 
ample, from NH 2 D observations we find clear evidence 
of the peak being offsets from the bright cores showing 
star formation activity (Pillai et al. 2011). The colder 
secondary cores on the other hand are expected to have 
abundant deuterated species. Are these secondary cores 
principally seen in H 2 D + emission representative of the 
prestellar/cluster gas? 

5.1. Extent o/H 2 D+ emission 

The most surprising element of our observation is the 
wide-spread emission of H 2 D + with no direct counter- 
part in previous dust continuum or line measurements 
(continuum data from Motte et al. 2007 and line data 
from Csengeri et al. 2011). The H 2 D + emission is in- 
stead distributed in cores around the bright continuum 
cores over a ~ 1 pc region. The individual cores have 
on average diameter of 0.2 pc roughly estimated by fit- 
ting ellipses by eye (34000 AU, P4 and P5 are smaller). 
This is more than twice the extent of H 2 D + emission ob- 
served towards the only two regions previously mapped 
in H 2 D + ; the low-mass prcstellar core in Taurus (L1544; 
Vastcl et al. 2006) and the low-mass cluster-forming core 
in Ophiucus (OphB2; Friesen et al. 2010). 

In the past, H 2 D + was thought to trace the innermost 
regions (few 1000 AU) of a dense core where all the heav- 
ier molecules deplete on to grain mantles with density 
exceeding 10 6 cm~ 3 (di Francesco et al. 2007). How- 
ever, with newer data on the collisional rate of H 2 D + 
the critical density of the transition has been corrected 
to ~ 10 5 cm~ 3 , an order of magnitude lower than previ- 
ously assumed (Caselli et al. 2008). Combined with the 



prevalent cold temperature in such dense cores, this im- 
plies that the H 2 D + distribution is expected to be much 
larger (» 7000 AU). The extent of H 2 D+ distribution in 
OphB2 or L1544 can then be comprehended. However, 
it is surprising that in Cyg-X this distribution is twice 
as large. Mapping of the entire DR21 filament is likely 
to discover more such cold quiescent gas pockets in an 
otherwise active high mass star-forming region. 

5.2. Absence of Dust Continuum peaks towards 
orf/io-H 2 D+ peaks 

Although there is weaker dust emission towards all 
H 2 D + cores, none of the H 2 D+ peaks have been detected 
as a peak in dust continuum emission with MAMBO as 
well as SMA. Depletion is often invoked to explain such 
lack of correlation between dust and gas in different trac- 
ers. However, the molecule in question is H 2 D + , one of 
the lightest deuterated ions that should survive deple- 
tion. Therefore, we suggest that dust continuum obser- 
vation is partially misleading, i.e. not tracing mass dis- 
tribution. We assert this by showing that H 2 D + emission 
comes from a region of low temperature and high density. 

We draw out the density in three different ways. First, 
let us define a spherical core around a H 2 D + core, for e.g. 
PI. It has a mass of 137 M within a radius of ~ 0.1 pc 
from MAMBO observations for a temperature of 15 K. 
This implies an average density of 5 x 10 6 cm~ 3 . Second, 
we extract the intensity of the filament in dust emis- 
sion towards an intensity minima in between two bright 
MAMBO cores, considerably offset from the H 2 D + cores 
and outside its 3 a boundary. Adopting the distance 
to the extreme ends of two H 2 D + cores boundary as a 
measure of the size of a typical core, in this case 60 arc- 
sec, we again derive density (2 x 10 5 cm -3 ). The crit- 
ical density of the observed H 2 D + transition at 20 K is 
~ 1.1 x 10 cm~ 3 . This shows that the density in a region 
larger than the H 2 D + core by a factor of ~ 60/20 = 3 
(20 arcsec corresponds to representative H 2 D + core size) 
is still high enough to get H 2 D + excited. Finally, N 2 H + 
4-3 emission is observed all along the filament including 
towards H 2 D + cores. The critical density 8 of this N 2 H + 
transition at 3 x 10 7 cm~ 3 is two orders of magnitude 
higher than that of the H 2 D+ transition, requiring very 
dense gas. Clearly not every region (like the envelope) 
in the filament can have such high densities or the im- 
plied high excitation temperature, and the line is likely 
sub-thermally excited (Evans 1999). We therefore used 
the non-LTE radiative transfer code RADEX 9 for a simple 
calculation of the expected N 2 H + 4-3 line brightness for 
a column density of 10 15 cm~ 2 and line width (FWHM) 
of 2 kins" 1 at 20 K for a wide range of densities (guided 
by N 2 H+ line parameters for the H 2 D+ cores given in Ta- 
ble 2-3). We are thus able to rule out densities below 
10 5 cm~ 3 for the H 2 D+ cores. 

Next, we look at the maximum temperature of the 
H 2 D + emitting region. H 2 D + cores are mid-IR-quiet 
at 24/xm, while the bright MAMBO cores have bright 
24/xm point like emission as observed with the MIPS 
instrument on Spitzer telescope (Spitzer Legacy Pro- 

8 Derived based on N2H+ collisional rate from LAMDA database 
for a temperature of 20 K (Schoier et al. 2005). 

9 van der Tak et al. (2007) 
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gram ID 40184, PI J. Hora), consistent with the notion 
that these bright objects are massive protostars (Chan- 
dler et al. 1993). A precise luminosity of these proto- 
stars is difficult to measure because the sources are satu- 
rated. Based on a single temperature Spectral Energy 
Distribution (SED) fit to their Far-Infrared(FIR) and 
1.3 mm data, Chandler et al. have determined the lumi- 
nosity of one of the protostars in the region (known as 
FIR1) to be ~ 1860 L© (this value has been corrected 
for our adopted distance of 1.7 kpc instead of 3 kpc). 
Given the poor spatial resolution of the IRAS data, the 
individual protostars would in any case not be distin- 
guished. We find two sources from the Red MSX Source 
(RMS) survey of a large sample of massive young stel- 
lar objects from the MSX point source catalog in our 
field (Hoare et al. 2005,Urquhart et al. 2008). The sur- 
vey database at www.ast.leeds.ac.uk/RMS characterises 
the luminosity of the two sources, G081. 7522+00.5906, 
G081. 7624+00.5916 as 2500 and 1000 L© respectively 
(Mottram et al. 2011). Assuming that FIR1/2/3 are in 
similar stages of evolution, we adopt 1000 L© as a con- 
servative estimate of the protostellar luminosity. The 
mass-weighted dust temperature within an aperture of 
radius R is approximated using Eq. (8) of Kauffmann 
et al. (2008). That equation assumes a simplified dust 
temperature dependence on radius and luminosity (i.e., 
calibrated by Terebey et al. 1993; T d oc L°- 2 /r 0A ), but 
further requires temperatures > 10 K. Assuming a den- 
sity profile g oc r -2 , these temperatures are then weighed 
by mass to obtain the mean temperature (e.g., Belloche 
et al. 2006). Following Eqs. (5), (8) from Kauffmann 
et al. (2008), the temperature at a radius of 30000 AU 
(~ 17 ") from the protostar would drop to 15 K. Here, 
the exact value of the radius is not relevant but within 
reasonable distance of the protostar the dust should be 
cold, a requisite for H2D + emission. 

So far we have shown that high density and low tem- 
perature — a necessary condition for H2D + excitation — 
is satisfied in an extended region around the protostars. 
Therefore, the cores seen in H 2 D+ need not necessarily 
correspond to a 3D maximum in density. To the extent 
of the validity of our first order calculations here in terms 
of effect or protostellar luminosity on heating, it is not 
required to have a density peak towards the H 2 D + peaks. 
If we were to do more elaborate calculations to derive the 
density structure, then a density peak would have en- 
hanced the cooling within the peak. However, the H2D+ 
peaks might well correspond to a temperature minima. 
We intend to investigate this via NH3 measurements. 

Why are the H2D+ cores not detected in SMA 1.3 mm 
dust continuum? The SMA continuum distribution is 
offset from the H2D 4 " peaks. Recall however that the 
H2D" 1 " peaks need not be density peaks. Therefore the 
fact that the dust peaks are offset is not inconsistent 
with the H2D+ observations. Even so, it is worth not- 
ing that all H 2 D + cores are widely associated with SMA 
dust continuum peaks, i.e. offset from rather than be- 
ing anti-correlated with continuum peaks. Based on the 
1.3 mm MAMBO data we calculated the total mass of 
the H2D" 1 " cores. We derive the gas mass at 15 K from 
the 1.3 mm following Kauffmann et al. (2008), using dust 
opacities at 1.3 mm (0.0102 cm 2 /g) for dust grains with 
thin ice mantles and gas density n(H) = 10 6 cm~ 3 as in 



Ossenkopf & Henning 1994, and gas to dust ratio of 100. 
The mass estimate is uncertain by a factor of 2 mainly 
due to the uncertain dust opacity. The typical mass of a 
H 2 D+ core is - 130 M . However, we find that < 10% 
mass resides in the cores detected with SMA. The miss- 
ing matter might be either dense but smooth structures 
filtered out by an interferometer like the SMA. This is 
validated by our N2D+ SMA data. It is contaminated by 
negative bowls (see third panel of Fig 1). This unwanted 
feature of synthesis array observations is a signature left 
by very extended objects which dominate the short spac- 
ings in the uv plane. Alternatively, a cluster of low mass 
cores well below the sensitivity limit of the array would 
also explain the missing matter. 

Given that 90% of the dust emission is missing in 
the SMA data, the lack of a one-to-one correlation be- 
tween SMA continuum and H 2 D + might not be signifi- 
cant. For example, consider the dust emission associated 
with clump PI. We know that ~ 90% of the dust emis- 
sion is filtered out with SMA. If the structure is smooth, 
then this ~ 90% emission could be randomly removed 
from any position in clump PI. Therefore the remain- 
ing emission might also be located at a random position 
within PI. Moreover, the offset between SMA dust and 
H2D + peaks is suggestive of density not being the only 
criteria to have H2D + peaks where they are observed. 
If it would have been the case, then naturally the peak 
with the highest density would be least likely to filtered 
out by SMA. Therefore we believe that the H2D 4 " peaks 
might be dominated by temperature effects. The dust 
might be significantly cooler than in the surroundings 
(and definitely cooler than at the bright dust continuum 
peaks, where protostars are embedded). A similar rea- 
soning was proposed by Di Francesco et al. (2004) who 
found that interferometric observations of N 2 H + actually 
betrayed the presence of a dense and cold starless core 
N6 in Oph A with no distinct peak observed in mm dust 
continuum observations. 

We have derived the virial mass of the H 2 D + cores 
following Bertoldi & McKee 1992. For a mean radius of 
0.1 pc and FWHM Av — 1.2 km s -1 , the virial mass of a 
H2D" 1 " core is 30 M Q . The average virial mass is 28 M Q , 
very small compared to the mass derived from MAMBO 
observations. Therefore these cores are gravitationally 
super-critical. The line width is highly supersonic (> 
1 kms -1 ), for all but one core (P4). The high line-width 
is another indication of a clumpy substructure within 
each H 2 D+ emitting region. 

5.3. Absence o/N 2 D+ peaks towards ort/io-H 2 D+ peaks 

As noted before, while N 2 D + emission as observed at 
high angular resolution with SMA is widely associated 
with H2D + , the peaks of these two ions however do not 
coincide. Lack of N2D + emission towards P2 is at least 
partially due to the problem of short spacing discussed in 
Section 5.2. Briefly, if the N 2 D + emission is as extended 
as the observed H2D + emission, then a significant frac- 
tion of this emission is likely to be filtered out by the 
SMA. We intend to pursue single dish observations to 
address this issue. Two cases deserve more attention 
and are discussed here. First, N2D 4 " is seen to skirt the 
H2D" 1 " emission in non-uniform clumpy structures except 
in the case of P2. The critical density for the 3-2 transi- 
tions of this line is an order of magnitude higher than our 
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H2D + transition. This implies that higher density is a 
requisite for its excitation. We have also shown that the 
H2D" 1 " peaks at low resolution need not be density peaks. 
The SMA continuum cores however should point us to- 
wards the highest density compact cores in the vicinity 
of H2D+ emission. It is exactly here that N2D + emission 
is concentrated. 

Next, two of the N 2 D+ clumps in the south-west part 
of the region are devoid of H 2 D + emission. It is dim- 
cult to understand why there are N 2 D + cores with no 
H 2 D + emission. Here, however excitation might play a 
role. Also, given that the SMA data is prone to arti- 
facts due to the uncertainties inherent in interfcrometric 
observations, these structures might be spurious. How- 
ever, we downloaded published and archived H 13 CO + 1- 
PdBI taken at similar resolution to our SMA observa- 
tions. The region is known as MDC N53 (Csengeri et al. 
2011). Clearly, there is H 13 CO + emission towards all the 
N2D + features including the two features not detected in 
H 2 D+. The peak N 2 D+ and H 13 CO+ line brightness are 
however much weaker (more than a factor of 4) toward 
these N 2 D+ cores. If the same scaling of line brightnesses 
is applied to H 2 D + , then it is evident why H 2 D+ is not 
detected (our highest H 2 D+ S/N is 8). 

5.4. Depletion of N-bearing molecules 

We will now discuss the distribution of N 2 H + and 
H 2 D + . At the low temperature of molecular clouds 
(T< 20 K), H 2 D + is the primary ion for deuterium 
transfer to other molecules including N2, which produces 
N2D + . Therefore, colder cores exhibit high deuteration 
as seen in the fractionation of molecules like N 2 H + , i.e 
high [N 2 D + ]/[N 2 H + ]. However, at such low tempera- 
tures molecules freeze out onto dust grains. This not 
only includes CO and its derivatives like HCO + but also 
N-bearing molecules like N2 (and consequently N 2 H + 
and NH3 , which are both formed by molecular nitro- 
gen). H2D + on the other hand can become abundant 
under such conditions where all heavy species have dis- 
appeared (depleted). To test this, we looked at the abun- 
dance distribution of N 2 H + over the region. While the 
average N 2 H + abundance of H 2 D + cores is 1 x 10~ 8 , 
no significant variation from those position offset from 
H2D + cores is found. The highest abundance is found 
towards the MAMBO bright peaks which are however 
higher by only a factor of ~ 2. Given that our resolu- 
tion is no finer than 34000 AU, this result is not surpris- 
ing. Stronger depletion might well be observed on much 
smaller spatial scales. While the estimated N2H + abun- 
dance of ~ 10 ~ 8 might seem significantly higher than 
that found in low-mass star forming regions with typical 
abundance of ~ 10~ 10 (Caselli et al. 2002a, Keto et al. 
2004) , we caution that the somewhat arbitrary choice of 
Tex (5 K) influences this calculation. The higher column 
densities characteristic of high mass star forming regions 
might explain the high abundance, however a Tex higher 
by 2 K is sufficient to lower the derived abundance by one 
order of magnitude. 

5.5. Deuterium Chemistry 

How does the observed ortho-H 2 D + abundance com- 
pare to existing chemical models? Given the poor spatial 
resolution of our H 2 D + observations, a tailored chemi- 
cal model is not warranted. Qualitative comparisons can 



still be made with models that have focused on the H 2 D + 
chemistry in dark clouds (Walmsley et al. 2004, Flower 
et al. 2004). To do this, we need to convert the observed 
ortho-H 2 D + abundance to the total H 2 D + abundance 
assuming an ortho-to-para (o/p) ratio for H 2 D + . Few 
observations have constrained this value in dark clouds 
(~ 0.8 in L1544 Vastcl et al. 2006). Assuming an o/p- 
H 2 D + ratio as in L1544, we get a total H 2 D+ abundance 
of ~ 3 x 10~ n . However, if the gas temperature is even a 
few degree warmer than in low-mass star forming regions, 
then the o/p-H 2 D+ would drop - as shown in Figure 6 of 
Flower et al. (2004). This is mainly due to dependence 
of o/p-H 2 D + on o/p-H 2 , and an increase in ortho-H 2 as 
shown by chemical models (Pagani et al. 1992). Such a 
drop in the ortho-H 2 D + column density was already no- 
ticed by Caselli et al. (2008) toward warmer starless cores 
and protostellar cores in their sample with higher deu- 
terium fractionation (Rn). We do not have temperature 
or i?D measurements towards all the H 2 D + cores. How- 
ever, Pillai et al. (2007) estimate a high rotational tem- 
perature of 18 K in the vicinity of P4 (within 40") and 
i?D = 0.34 where Rr> is measured as the ratio of NH 2 D 
to NH 3 column density. This corresponds to a o/p-H 2 D+ 
ratio of about 0.6 (from Fig. 6 of Flower et al. 2004). 
We also compared our ortho-H 2 D + abundance to the dy- 
namical model of Flower et al. (2006) that incorporates 
deuterium chemistry in the early stages of protostellar 
collapse, specifically their Figure 6 where ortho-H 2 D+ 
abundance at 10 K is plotted as a function of the gas 
density during collapse. For a density of ~ 5 x 10 5 cm~ 3 
o/p-H 2 close to 0.03 is consistent with our observations. 
An o/p-H 2 of << 1 that we indirectly infer is consistent 
with those measured in nearby low mass star forming 
regions. Observations of deuterated molecules with its 
parent molecule for e.g. DCO+/HCO+ (Maret & Bergin 
2007), N 2 D + N 2 H + (Pagani et al. 2009), non-LTE radia- 
tive transfer modeling of H 2 CO absorption (Troscompt 
et al. 2009), abundance ratio of Nitrogen hydrides (Dis- 
laire et al. 2012) have all delivered o/p-H 2 of << 0.1. 
This implies that the para-H 2 dominates in molecular 
clouds. Similar work in high mass star forming region is 
yet to be done. Detailed modeling should be made with 
higher resolution data together with tighter constraints 
on the density and temperature in Cygnus-X. 

5.6. H 2 D + emission in high and low-mass star forming 

regions 

We present here the first map of H 2 D + distribution 
in a high mass star forming region. Previously, single 
spectra of H2D + have been reported in at least two such 
regions, OriB9, a 100 M© cold, quiescent core in Ori- 
onB (Harju et al. 2006) and ~ 4 a detection in a core 
in IRDC MSX G030.88+00.13 (Swift 2009). The ortho- 
H 2 D + abundance in Cygnus-X is in between the values 
for OriB9 (~ 10~ 10 ) and IRDCs - 10~ 13 . However, 
Cygnus-X is also intermediate between these two regions 
in terms of its distance from us. Therefore, the extreme 
range in H 2 D+ abundance might merely reflect the ef- 
fect of beam dilution. However, the line width of our 
H2D" 1 " cores are similar to the core in G030. 88+00. 13 
(0.9 kms -1 ), and a factor of two higher than in OriB9 
(~ 0.4 kms -1 ). We also note that due to a lack of con- 
tinuum data, Harju et al. derive OriB9 abundance based 
on the known Ammonia abundance in the region (not to- 
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wards the core), and might be offset from its true value. 
The H2D + column density in OrionB9 itself is compa- 
rable to Cygnus-X. The difference in line width is also 
reconciled if we account for the different spatial scales 
probed. Using the line width-size relation with an expo- 
nent of 0.4, if we were to scale the line width at a physical 
resolution of 0.03 pc i.e. 17 " FWHM APEX beam at 
410 pc to that at 0.16 pc (13 " FWHM JCMT beam at 
1700 pc), the expected line width is even higher than the 
observed line width. Therefore the physical conditions 
in Cygnus-X might be similar to those in the quiescent 
and massive core in Orion. 

As noted before, the average ortho-H2D + abundance 
in Cygnus-X is a factor of 1- 20 lower than that in low- 
mass cores (Caselli et al. 2008). Assuming that the true 
abundance is similar to what is measured in low-mass 
star forming regions, this might be due to a low filling 
factor. In that case H 2 D+ emission would be widespread 
but concentrated in small regions, with a filling factor 
of ~ 0.1. H2D + has been mapped in two low- mass star 
forming regions, Oph B2, a dense core in the low-mass 
cluster forming region in Ophiucus (Friesen et al. 2010), 
and L1544, a prestellar core in Taurus (Vastel et al. 
2006). The peak H2D + abundance in L1544 is an or- 
der of magnitude higher than for Cygnus-X. Vastel et 
al. have used a spherically symmetric chemical model 
incorporating the physical structure of the cloud (den- 
sity, temperature), and fit the abundance variation of 
H 2 D + as a function of distance from the core center. 
While the best fit to their H 2 D+ finds a peak at lO" 10 , 
Xo-h 2 d+ drops to ~ 10~ n at ~ 10000 AU from the cen- 
ter, still smaller than our typical ~ 17000 AU radii. This 
re-affirms our earlier conclusion that these observations 
have uncovered extended cold dense gas that were un- 
detected with typical molecular line and dust continuum 
surveys. Motte et al. (2007) surveyed the entire complex 
with MAMBO with the goal of identifying high mass 
cores at the earliest stage of evolution, but did not find 
high-mass analogs of prestellar dense cores. In the region 
of overlap with our H 2 D+ map, a single object (CygX- 
N53) has been identified as a massive dense core, which 
is distinctly offset from the H2D+ cores. However, the 
three massive dense, colder (and younger) cores discov- 
ered by this work (PI, P2, P3) are also undetected in line 
observations (H 13 CO+ and HCN 1-0) of Csengeri et al. 
2011. 

Now, L1544 is a prestellar core with just enough mass 
(8 M Q , Tafalla et al. 1998) to form a single low-mass star. 
Oph B2 on the other hand is a low-mass cluster forming 
core. Many of the H2D + line properties in Cygnus-X are 
mirrored in Oph B2. Namely, the H 2 D + peak is offset 
from the dust as well as N 2 D + peaks. The line widths are 
higher than in L1544, although still transonic. The emis- 
sion is also distributed over a larger region. H2D + core 
sizes are not constrained by these observations, however 
the whole emission extends to ~ 14000 AU. Naively, the 
more extreme H2D + properties in a high mass region like 
Cygnus-X seems to be scaled up from low-mass clusters. 
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6. SUMMARY 

We present the first map of ort/io-H 2 D+ J fe + fc - = 1^ 
— > li ; i transition in a high mass star forming region. 
Distributed over a physical scale of 1 pc in at least five 
clumpy structures, this is also the largest H2D + map of 
any star forming region. The detection is weak and shows 
an anti-correlation with star-formation, wherein H2D+ 
clearly avoids the active star forming clusters, however is 
still in their vicinity and has no IR emission of its own. 
The H2D" 1 " emission is widely associated with dense dust 
emission as well as N2D + as observed with SMA at high 
angular resolution. However, surprisingly the dust and 
N 2 D + peaks are distinctly offset from the H 2 D+ peaks. 
While the derived H 2 D + abundance is up to an order of 
magnitude lower than in low-mass star forming regions, 
the N 2 D+ and N 2 H+ column density and abundance is 
comparable to or even higher than low-mass star forming 
regions. The lack of local dust counterparts remains to 
be better understood. Our best reasoning is that H 2 D+ 
cores might have significantly clumpy sub-structure with 
temperature minima at the observed H2D 4 " peaks that 
need not be density peaks. We have thus shown that 
H 2 D + is as good and unique a tool in tracing the cold 
and dense gas in low-mass star forming regions as in such 
complex massive star forming regions. We have discov- 
ered a large reservoir of cold and dense gas undetected 
by previous surveys that were targeted at identifying the 
earliest stages of high-mass star formation. This suggests 
that submm continuum studies alone which are aimed at 
detecting massive prestellar cores might fail to identify a 
large population of such cores in earliest stages. Sensi- 
tive H2D" 1 " mapping of the entire DR21 filament is likely 
to discover more of such cold quiescent gas reservoirs in 
an otherwise active high mass star-forming region. 
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Fig. 1. — Panel A (from left): SMA 1.3mm dust continuum map which is a mosaic of four pointings (color scale), and MAMBO 1.3 mm 
map (contours, Motte et al. 2007). Panel B: MAMBO 1.3 mm dust continuum and H2D+ S/N map as contours (contours, 3er in steps of 
1.5er where 0.013 < cr/K < 0.021). Panel C: N2D+ 3-2 moment map (weighted mean intensity integrated from -14 kms^ 1 to -8 kms -1 .), 
and JCMT H2D+ S/N map (contours, 3cr in steps of 1.5<r). Panel D: N2H+ 4-3 moment map (weighted mean intensity integrated from 
-10 kms -1 to 3 kms -1 , grey scale), and JCMT H2D+ S/N map (contours, 3cr in steps of 1.5cr). The H2D+ core positions are marked as 
P'n'. Stars denote the embedded population identified by Kumar et al. (2007). 
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Fig. 2.— Left Panel: JCMT spectra of H 2 D+ at 372.421 GHz. Centre Panel: JCMT N 2 H+ 4-3 at 372.672 GHz. Right Panel: SMA 
N 2 D+ 3-2 at 231.5 GHz. Spectra for all tracers have been extracted for the same spatial resolution, and the offset positions (in arcsecs) are 
as indicated in the figure. Note that the SMA spectra have a factor of 2 poorer spectral resolution than the JCMT spectra (~ 0.2 kms -1 ). 
H 2 D+ LSR velocity for each position is indicated as the line of reference in the three panels. 



TABLE 2 

Gaussian fit parameters from JCMT Observations. 
Source (RA, Dec.) Molecule Tmb Vlsr At) Ai> nt Velocity Range j Tmb dv 



J2000 




K kms" 1 


kms 1 


kms 


kms 1 


K kms" 1 


PI H 2 D+ 


ll,0" 


-li,i0.18 


(0.02)-4.6 


(O.i; 


1.20 


(0.08) 


1.13 






(20:39:03.7,+42:25:21)N 2 H+ 


4-3 


3.27 


(0.03)-4.2 


(o.r 


1.77 


(0.02) 


1.76 


(-8.0,0.0) 


9.65 (0.06) 


N 2 D+ 


3-2 


0.73 


(0.14)-4.8 


(0.1 


0.82 


(0.14) 


0.81 






P2 H 2 D+ 


li,o- 


-li,i0.14 


(0.02)-4.6 


(0.1 


1.50 


(0.17) 


1.44 






(20:39:03.1,+42:25:37)N 2 H+ 


4-3 


4.15 


(0.02)-4.2 


(0.1 


2.04 


(0.01) 


2.03 


(-8.0,0.0) 


13.72 (0.03) 


N 2 D+ 


3-2 


0.24 


(0.12)-4.7 


(0.2; 


0.85 


(0.49) 


0.84 






P3 H 2 D+ 


li,o- 


-li,i0.15 


(0.02)-4.3 


(0.1; 


1.26 


(0.13) 


1.19 






(20:39:01.7, +42:25:21) N 2 H+ 


4-3 


3.69 


(0.02)-3.7 


(0.1; 


1.69 


(0.16) 


1.68 


(-8.0,0.0) 


10.79 (0.04) 


N 2 D+ 


3-2 


0.33 


(0.13)-3.5 


(0.1; 


0.67 


(0.23) 


0.65 






P4 H 2 D+ 


li.o- 


-lnO.12 


(0.02)-2.8 


(0.1; 


1.61 


(0.15) 


1.56 






(20:38:58. 3, +42:24:44) N 2 H+ 


4-3 


2.40 


(0.02)-2.6 


(0.1; 


1.47 


(0.01) 


1.46 


(-8.0,2.0) 


6.93 (0.03) 


N 2 D+ 


3-2 


0.19 


(0.07)-2.2 


(0.1; 


0.69 


(0.99) 


0.67 






P5 H 2 D+ 


li.o- 


-li,i0.13 


(0.07)-3.6 


(0.1; 


0.82 


(0.35) 


0.71 






(20:38:58.9+42:24:22) N 2 H+ 


4-3 


2.77 


(0.02)-3.4 


(0.1 


1.92 


(0.01) 


1.91 


(-8.0,2.0) 


9.30 (0.04) 


N 2 D+ 


3-2 


0.20 


(0.08)-3.3 


(0.1; 


0.74 


(0.28) 


0.72 







Av and Ati n t are the observed and non-thermal FWHM line width, respectively. Av n t is calculated after accounting for thermal line width 
at 15 K. Towards position P4, the H 2 D+ spectrum shows an asymmetry which if taken as due to two velocity components has line width 
(FWHM) of 0.88 and 0.68 kms _1 at LSR velocity of -3.3 and -2.4 kms -1 , respectively. The FWHM derived for H 2 D+ is from a simple 
Gaussian fit while that for N 2 D" 1 " and N 2 H + is from a hyperfine fit. 



H2D+ in a high mass star- forming region 
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TABLE 3 

H2D+, N2H+COLUMN DENSITY AND ABUNDANCE. 



Position 


N(o - H 2 D+) a 


X -h 2 d+ b 


N(N 2 H+) C 


XN 2 H+ d 




(xlO 12 ) cm" 2 


(xlO" 11 ) 


(xlO 15 ) cm" 2 


(XlO" 8 ) 


PI 


3.8 (0.5) 


1.9 


2.0 (0.01) 


1.0 


P2 


3.7 (0.7) 


1.4 


2.8 (0.01) 


1.1 


P3 


3.3 (0.6) 


1.7 


2.2 (0.01) 


1.1 


P4 


3.4 (0.7) 


2.1 


1.4 (0.01) 


0.9 


P5 


1.9 (1.3) 


0.9 


1.9 (0.01) 


0.9 



a The column density determined for excitation temperature, Tex= 10 K. 

b Fractional orffto-H2D+abundance based on H2 column density derived from MAMBO 1.3 mm data for dust temperature, Tkin= 15 K. 
c The N2H + column density determined for excitation temperatures, Tgx= 5 K. 

d Fractional N2H+abundance based on H2 column density derived from MAMBO 1.3 mm data for dust temperature, Tkin= 15 K. 



